function FilterNature(n2)

for n1=1860:2021
    eval(sprintf('VecResultsArticle%d={};',n1));
    eval(sprintf('VecResultsArticleType%d={};',n1));
end
ValueFolder='/sps/isc/rneveu/EtudeIRM/Nature/Results/';
ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
VecArticleType={'Academic Dissertation' 'Adaptive Clinical Trial' 'Clinical Trial' 'Comparative Study' 'Controlled Clinical Trial' 'Equivalence Trial' 'Evaluation Study' 'Journal Article' 'Meta-Analysis' 'Multicenter Study' 'Observational Study' 'Pragmatic Clinical Trial' 'Randomized Controlled Trial' 'Review' 'Statistics' 'Technical Report' 'Systematic Review' 'Twin Study' 'Validation Study'};
VecArticleType2={'Review' 'Systematic Review'};
VecA=192:1:198;
VecA=[VecA 65];
Veca=224:1:230;
Veca=[Veca 97];
VecC=[199 67];
Vecc=[231 99];
VecE=200:203;
VecE=[VecE 69];
Vece=232:1:235;
Vece=[Vece 101];
VecI=204:207;
VecI=[VecI 73];
Veci=236:1:239;
Veci=[Veci 105];
Vecn=[241 110];
VecO=210:214;
VecO=[VecO 216 79];
Veco=242:1:246;
Veco=[Veco 248 111];
VecU=217:1:220;
VecU=[VecU 85];
Vecu=249:1:252;
Vecu=[Vecu 117];
VecY=[221 89];
Vecy=[253 121];
Vecs=[223 115];
VecSpecialCharacters=[{VecA} {Veca} {VecC} {Vecc} {VecE} {Vece} {VecI} {Veci} {Vecn} {VecO} {Veco} {VecU} {Vecu} {VecY} {Vecy} {Vecs}];
ValueClock=clock;
disp(['Start of the processing at ' num2str(ValueClock(1,4)) ':' num2str(ValueClock(1,5)) ':' num2str(ValueClock(1,6)) ' for journal ' num2str(n2)]);
ValueFile=['ResultsJournal' num2str(n2) '.mat'];
ValueFileTemp=dir([ValueFolderTemp ValueFile]);
if size(ValueFileTemp,1)==0
    n1Temp2=0;
    try
        ValueFile=['ResultsJournalNatureTemp' num2str(n2) '.mat'];
        load([ValueFolderTemp ValueFile]);
        n1Temp2=n1;
    end
    n1Temp2=n1Temp2+1;
    for n1=n1Temp2:1062
        n1
        eval(sprintf('load([ValueFolder ''ResultsFile%d.mat'']);',n1));
        eval(sprintf('VecValueTemp=VecResult2Journal%d;',n2));
        for n3=1:size(VecValueTemp,1)
            ValueTemp=VecValueTemp{n3,1};
            if size(ValueTemp,2)>1
                [ValueFirst,ValueLast,ValueAffiliation,ValueORCID,nAuthors]=ExtractAuthor(ValueTemp,0,1);
                ValueArticleType=ComputeArticleType(VecValueTemp{n3,9},VecArticleType,VecValueTemp{n3,10},VecValueTemp{n3,11});
                ValueArticleType2=ComputeArticleType(VecValueTemp{n3,9},VecArticleType2,VecValueTemp{n3,10},VecValueTemp{n3,11});
                if ValueArticleType>0
                    if nAuthors>=1
                        ValueArticleType3=1;
                    else
                        ValueArticleType3=2;
                    end
                else
                    if ValueArticleType2>0
                        ValueArticleType3=3;
                    else
                        ValueArticleType3=4;
                    end
                end
                ValueDate=VecValueTemp{n3,5};
                [ValueYear,ValueMonth,ValueDay]=ExtractDate(ValueDate);
                if size(ValueYear,2)>0
                    ValueDateArticleTemp=[{ValueYear} {ValueMonth} {ValueDay}];
                    ValueYearTemp=str2num(ValueYear);
                    eval(sprintf('VecResultsArticle=VecResultsArticle%d;',ValueYearTemp));
                    eval(sprintf('VecResultsArticleType=VecResultsArticleType%d;',ValueYearTemp));
                    VecResultsArticleType(size(VecResultsArticleType,1)+1,:)=[{size(VecResultsArticleType,1)+1} {nAuthors+1} {ValueArticleType} {ValueArticleType2} {ValueArticleType3} {nAuthors+1} {ValueYear} {ValueMonth} {ValueDay}];
                    n3Temp=size(VecResultsArticle,1);
                    VecResultsArticleTemp=[];
                    VecAffiliationTemp={};
                    ValueAffiliationTemp='';
                    nAffiliation=0;
                    VecResultsTotalTemp=zeros(1,6);
                    VecResultsOriginalTemp=zeros(1,6);
                    VecResultsReviewTemp=zeros(1,6);
                    for n41=1:nAuthors
                        [ValueFirst,ValueLast,ValueAffiliation,ValueORCID,nAuthorsTemp]=ExtractAuthor(ValueTemp,n41,2);
                        VecAffiliationTemp(1,n41)={ValueAffiliation};
                        if ValueArticleType>0&&nAuthors>1
                            VecResultsOriginalTemp(1,1)=1;
                            VecResultsTotalTemp(1,1)=1;
                            VecResultsOriginalTemp(1,6)=n41;
                            VecResultsTotalTemp(1,6)=n41;
                            if n41==1
                                VecResultsOriginalTemp(1,2)=1;
                                VecResultsTotalTemp(1,2)=1;
                            else
                                if n41<nAuthors
                                    if n41==2
                                        VecResultsOriginalTemp(1,3)=1;
                                        VecResultsTotalTemp(1,3)=1;
                                    else
                                        VecResultsOriginalTemp(1,4)=1;
                                        VecResultsTotalTemp(1,4)=1;
                                    end
                                else
                                    VecResultsOriginalTemp(1,5)=1;
                                    VecResultsTotalTemp(1,5)=1;
                                end
                            end
                        else
                            if ValueArticleType2>0&&nAuthors>1
                                VecResultsReviewTemp(1,1)=1;
                                VecResultsTotalTemp(1,1)=1;
                                VecResultsReviewTemp(1,6)=n41;
                                VecResultsTotalTemp(1,6)=n41;
                                if n41==1
                                    VecResultsReviewTemp(1,2)=1;
                                    VecResultsTotalTemp(1,2)=1;
                                else
                                    if n41<nAuthors
                                        if n41==2
                                            VecResultsReviewTemp(1,3)=1;
                                            VecResultsTotalTemp(1,3)=1;
                                        else
                                            VecResultsReviewTemp(1,4)=1;
                                            VecResultsTotalTemp(1,4)=1;
                                        end
                                    else
                                        VecResultsReviewTemp(1,5)=1;
                                        VecResultsTotalTemp(1,5)=1;
                                    end
                                end
                            else
                                VecResultsTotalTemp(1,1)=1;
                                VecResultsTotalTemp(1,6)=n41;
                                if n41==1
                                    VecResultsTotalTemp(1,2)=1;
                                else
                                    if n41<nAuthors
                                        if n41==2
                                            VecResultsTotalTemp(1,3)=1;
                                        else
                                            VecResultsTotalTemp(1,4)=1;
                                        end
                                    else
                                        VecResultsTotalTemp(1,5)=1;
                                    end
                                end
                            end
                        end
                        if size(ValueAffiliation,2)>1&&size(ValueAffiliationTemp,2)==0
                            ValueAffiliationTemp=ValueAffiliation;
                            nAffiliation=nAffiliation+1;
                        else
                            if size(ValueAffiliation,2)>1&&size(ValueAffiliationTemp,2)>1
                                ValueResult=CompareAffiliation(ValueAffiliation,ValueAffiliationTemp);
                                if ValueResult==0
                                    nAffiliation=nAffiliation+1;
                                end
                            end
                        end
                        TestTemp3=0;
                        n43=1;
                        while n43<=n3Temp&&TestTemp3==0
                            if strcmps(double(VecResultsArticle{n43,1}),double([ValueLast ' ' ValueFirst]),VecSpecialCharacters)>=0.99
                                TestTemp3=n43;
                            end
                            n43=n43+1;
                        end
                        if TestTemp3==0
                            VecResultsArticle(n3Temp+1,1)={[ValueLast ' ' ValueFirst]};
                            VecResultsArticle(n3Temp+1,2)={[n1 n3]};
                            VecResultsArticle(n3Temp+1,3)={VecResultsTotalTemp};
                            VecResultsArticle(n3Temp+1,4)={VecResultsOriginalTemp};
                            VecResultsArticle(n3Temp+1,5)={VecResultsReviewTemp};
                            VecResultsArticle(n3Temp+1,6)={{ValueAffiliation}};
                            VecResultsArticle(n3Temp+1,7)={{ValueDateArticleTemp}};
                            VecResultsArticleTemp=[VecResultsArticleTemp n3Temp+1];
                        else
                            ValueTemp4=VecResultsArticle{TestTemp3,2};
                            ValueTemp4=[ValueTemp4;n1 n3];
                            VecResultsArticle(TestTemp3,2)={ValueTemp4};
                            ValueTemp4=VecResultsArticle{TestTemp3,3};
                            ValueTemp4=[ValueTemp4;VecResultsTotalTemp];
                            VecResultsArticle(TestTemp3,3)={ValueTemp4};
                            ValueTemp4=VecResultsArticle{TestTemp3,4};
                            ValueTemp4=[ValueTemp4;VecResultsOriginalTemp];
                            VecResultsArticle(TestTemp3,4)={ValueTemp4};
                            ValueTemp4=VecResultsArticle{TestTemp3,5};
                            ValueTemp4=[ValueTemp4;VecResultsReviewTemp];
                            VecResultsArticle(TestTemp3,5)={ValueTemp4};
                            ValueTemp4=VecResultsArticle{TestTemp3,6};
                            ValueTemp4=[ValueTemp4;{ValueAffiliation}];
                            VecResultsArticle(TestTemp3,6)={ValueTemp4};
                            ValueTemp4=VecResultsArticle{TestTemp3,7};
                            ValueTemp4=[ValueTemp4;{ValueDateArticleTemp}];
                            VecResultsArticle(TestTemp3,7)={ValueTemp4};
                            VecResultsArticleTemp=[VecResultsArticleTemp TestTemp3];
                        end
                    end
                    if nAffiliation==1
                        for n44=1:size(VecAffiliationTemp,2)
                            if size(VecAffiliationTemp{1,n44},2)>=1
                                ValueAffiliation=VecAffiliationTemp{1,n44};
                            end
                        end
                        for n44=1:size(VecAffiliationTemp,2)
                            VecAffiliationTemp(1,n44)={ValueAffiliation};
                        end
                    end
                    for n44=1:size(VecAffiliationTemp,2)
                        ValueTemp4=VecResultsArticle{VecResultsArticleTemp(1,n44),6};
                        ValueTemp4(end,1)=VecAffiliationTemp(1,n44);
                        VecResultsArticle(VecResultsArticleTemp(1,n44),6)={ValueTemp4};
                    end
                    eval(sprintf('VecResultsArticle%d=VecResultsArticle;',ValueYearTemp));
                    eval(sprintf('VecResultsArticleType%d=VecResultsArticleType;',ValueYearTemp));
                end
            end
        end
        if mod(n1,50)==0
            ValueFile=['ResultsJournalNatureTemp' num2str(n2) '.mat'];
            save([ValueFolderTemp ValueFile]);
        end
    end
    ValueTemp=['save([ValueFolderTemp ''ResultsJournal' num2str(n2) '.mat''],'];
    for n1=1860:2021
        ValueTemp=[ValueTemp '''VecResultsArticle' num2str(n1) ''',' '''VecResultsArticleType' num2str(n1) ''','];
    end
    ValueTemp(1,end)=')';
    ValueTemp=[ValueTemp ';'];
    eval(sprintf(ValueTemp));
end
ValueClock=clock;
disp(['End of the processing at ' num2str(ValueClock(1,4)) ':' num2str(ValueClock(1,5)) ':' num2str(ValueClock(1,6)) ' for journal ' num2str(n2)]);